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Abstract. We review the approach to modeling soft hadron physics observables based on the Dyson-Schwinger equations 
of QCD. The focus is on light quark mesons and in particular the pseudoscalar and vector ground states, their decays and 
electromagnetic couplings. We detail the wide variety of observables that can be correlated by a ladder-rainbow kernel with one 
infrared parameter fixed to the chiral quark condensate. A recently proposed novel perspective in which the quark condensate 
is contained within hadrons and not the vacuum is mentioned. The valence quark parton distributions, in the pion and kaon, 
as measured in the Drell Yan process, are investigated with the same ladder-rainbow truncation of the Dyson-Schwinger and 
Bethe-Salpeter equations. 
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DYSON-SCHWINGER EQUATIONS OF QCD 

A great deal of progress in the QCD modeling of hadron physics has been achieved through the use of the ladder- 
rainbow truncation of the Dyson-Schwinger equations (DSEs). The DSEs are the equations of motion of a quantum 
field theory. They form an infinite hierarchy of coupled integral equations for the Green's functions (n-point functions) 
of the theory. Bound states (mesons, baryons) appear as poles in the appropriate Green's functions, and, e.g., the Bethe- 
Salpeter bound state equation appears after taking residues in the DSE for the appropriate color singlet vertex. The 
Faddeev equation for the 3-quark structure of baryons follows in a similar way from the residues of the appropriate 
vertex DSE. For recent reviews on the DSEs and their use in hadron physics, see Refs. [1, 2, 3, 4]. 
In the Euclidean metric that we use throughout, the DSE for the dressed quark propagator is 

S{p)- 1 = Z2ip , +Z A m{pL)+Z x J^g 2 D llv (p-q)—Y il S{q)T l v (q,p), (1) 

where D^ v (k) is the renormalized dressed-gluon propagator, F' v (q, p) is the renormalized dressed quark-gluon vertex. 
We use f£ to denote J h d A q/(2n) A with A being the mass scale for translationally invariant regularization. The 
renormalization condition is S(p)~ l = iy- p + m(jj.) at a sufficiently large spacelike jj, 2 , with m(fl) the renormalized 
mass at renormalization scale jj,. We use jj, = 19GeV. The Z,(jU, A) are renormalization constants. 

Bound state pole residues of the inhomogeneous Bethe-Salpeter equation (BSE) for the relevant vertex, yield the 
homogeneous BSE bound state equation 

r ai ( P+ ,P-) = [ A K(p,q;P)S a (q + )r ah (q + ,q_)S h (q-) , (2) 
Jq 

where K is the renormalized qq scattering kernel that is irreducible with respect to a pair of qq lines. Quark momenta 
are q + = q + rjP and q_ = q— (1 — Tj)P where the choice of 77 is equivalent to a definition of relative momentum q; 
observables should not depend on 77. The meson momentum satisfies P 2 = —M 2 . 

A viable truncation of the infinite set of DSEs should respect relevant (global) symmetries of QCD such as 
chiral symmetry, Lorentz invariance, and renormalization group invariance. For electromagnetic interactions and 
Goldstone bosons we also need to respect color singlet vector and axial vector current conservation. The rainbow- 
ladder (LR) truncation achieves these ends by the replacement K(p,q;P) — > — 4n a e ff(k 2 ) D f ^ (k) ^ 7^ ® ^Yv along 
with the replacement of the DSE kernel by Z 1 g 2 D llv (k)F' v (q, p) A%a e {f{k 2 )D i ^{k)y v ^- where k = p - q, and 

<*eff(£ 2 ) is an effective running coupling. This truncation is the first term in a systematic expansion [5, 6] of K; 
asymptotically, it reduces to leading-order perturbation theory. These two truncations are mutually consistent: together 



TABLE 1. DSE results [10, 11, 12, 13, 14, 15] for pseudoscalar and vector meson masses and 
electroweak decay constants, together with experimental data [16]. Units are GeV except where 
indicated. Quantities marked by f are fitted with the indicated current quark masses and the infrared 
strength parameter of the ladder-rainbow kernel. 
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they produce color singlet vector and axial-vector vertices satisfying their respective Ward identities. This ensures that 
the chiral limit ground state pseudoscalar bound states are the massless Goldstone bosons from dynamical chiral 
symmetry breaking [7, 8]; and ensures electromagnetic current conservation [9]. The ladder-rainbow kernel found to 
be successful for light quark hadrons [8, 10] can be written a e ff(£ 2 ) = a 1R (k 2 ) + a vv (k 2 ). The IR term implements 
the strong infrared enhancement in the region < k 2 < 1 GeV 2 required for sufficient dynamical chiral symmetry 
breaking. The UV term preserves the one-loop renormalization group behavior of QCD: a e ff(£ 2 ) — > (X s (k 2 ) n °° v (k 2 ) 
in the ultraviolet with Nf = 4 and Aqcd = 0.234GeV. The strength of a 1R along with two quark masses are fitted to 
(qq), m njK . 

Selected light quark meson results are displayed in Table 1 . The charge radii are obtained from the calculated charge 
form factors of % and K [11, 12] in the triangle diagram approximation. This approximation is the dynamical partner 
of LR truncation: charge is automatically conserved. Vector meson dominance is a natural outcome here since it is 
generated by the inhomogeneous BSE for the dressed photon-quark vertex. Related to this, the decay constants for 
a vector meson to a pair of pseudoscalars are obtained by direct computation of the triangle diagram [13], or from 
analysis of the timelike electromagnetic form factors [14]. The Abelian axial anomaly coupling g ny y and associated 
interaction radius are obtained from electromagnetic transition form factors [14], which also yield the radiative decay 
constants. The three-body leptonic decays of the kaon [15] are also shown. 

It has remained true for many years that efforts [ 1 7] to solve truncated versions of the DSEs for the QCD gauge sector 
n-point functions that determine the ladder-rainbow kernel have been unable to produce sufficient kernel strength. 
That is, an empirically acceptable amount of dynamical chiral symmetry breaking (and associated ground state 
hadron mass scales) has long been underestimated that way. Thus our infrared kernel component has remained semi- 
phenomenological while we seek to characterize it's deficiencies through application to new classes of observable. 
Some efforts to go beyond the semi-phenomenological ladder-rainbow truncation so as to progressively reduce the 
amount of phenomenology have employed an explicit 1-loop dressing of the gluon vertex: a gluon loop [18], or a 
pion loop [19]. However, an amount of phenomenology is still needed in these approaches to reproduce light quark 
physics. Furthermore, a study [20, 21] of this type of a diagram-by-diagram approach to the vertex and BSE kernel, 
within a schematic model of simplified momentum dependence, shows that the produced color factors run rapidly out 
of control. 

It is only recently that large volume lattice-QCD calculations [22] and analysis have established that the gluon 2- 
point function is non-zero in the IR limit, corresponding to a gluon effective mass of <~ 0.5 GeV, and that the ghost 2- 
point function is no more singular than its perturbative behavior. These modern lattice-QCD results have recently [23] 
been combined to produce a kernel for the quark DSE that yields an empirically acceptable amount of dynamical 
chiral symmetry breaking. The gluon-quark vertex in that approach might be used in the recently developed method 



for generalized BSE kernel construction [24] to produce a physical and capable pathway to a wider variety of hadron 
bound state calculations. 



NOVEL PERSPECTIVE ON THE QUARK CONDENSATE 

Recently it has been advocated that the quark condensate, rather than being a spacetime independent constant filling all 
space, is really a property of the Bethe-Salpeter wavefunction of hadrons [25]. Clear evidence for such a perspective 
is provided by the exact QCD expression for the residue p K at the pseudoscalar meson pole in the 3-point function 
describing the coupling of a pseudoscalar field with quarks [7] 

ip K = -{Q\qij 5 q\n) =Z 4 (C,A) tr cd f y 5 S(q + )T a (q;P)S(q^) . (3) 

Due to Goldsone's Theorem, and the resulting Goldberger-Treiman relation in which the chiral meson's invariant BSE 
amplitude associated with 75 is the chiral quark scalar self-energy divided by the chiral f n , one can show rigorously 
that [7] 

km f n (Q\qy 5 q\n) = -Z 4 (C,A)tr cd / 7^4 S Q (q;Q = (qq)\. (4) 

Thus the chiral condensate is a property of the qq projection of the chiral pion Bethe-Salpeter wavefunction, and thus 
from this perspective is not defined away from this finite-sized bound state. Although it is true that the (point) pion 
cloud effect on the pion leads to an infinite charge radius, the qq dressed quark core remains of finite size and that is 
where Eq. (4) places the condensate. Furthermore this is in accord with quark and gluon confinement; these virtual 
fields only exist inside their hadronic containers. 

The quark condensate represents an enormous energy density, comparable to the interior of compact stellar objects. 
The above perspective would eliminate the 10 45 overestimate of the cosmological constant from the QCD vacuum 
energy density. If the electroweak quanta turn out to be composites of higher mass entities, like technicolor quarks, the 
vacuum field we presently call the Higgs scalar is needed only inside those containers. This could eliminate the 10 56 
discrepancy from that sector of the Standard Model. 



VALENCE QUARK DISTRIBUTIONS IN THE PION AND KAON 

Data for the momentum-fraction probability distributions of quarks and gluons in the pion have primarily been inferred 
from Drell-Yan [26, 27, 28] and direct photon production [29] in pion-nucleon and pion-nucleus collisions, and semi- 
inclusive ep — > eNX reactions [30]. For a recent review of nucleon and pion parton distributions see Ref. [31]. 

Lattice-QCD is restricted to low moments of the distributions, not the distributions themselves [37]. Model calcu- 
lations of deep inelastic scattering (DIS) parton distribution functions are challenging because it is necessary to have 
perturbative QCD features (including the evolution of scale) coexisting with a covariant nonperturbative model made 
necessary by the bound state nature of the target. Aspects of chiral symmetry have led to DIS calculations within the 
Nambu-Jona-Lasinio model [38, 39, 40] but there are a number of difficulties with this approach [39, 40], among 
them a point structure for the pion BS amplitude at low model scale and a marked sensitivity to the regularization 
procedure due to the lack of renormalizability. Constituent quark models have also been employed [41, 42], with the 
difficulties encountered in such studies considered in Ref. [43]. Instanton-liquid models [44] have also been used. In 
these approaches, it is difficult to have pQCD elements join smoothly with nonperturbative aspects. All of these issues 
can in principle be addressed if parton distribution functions can be obtained from a model based on the DSEs, and in 
particular styled after the approach we have previously made for the related pion electromagnetic form factor [11, 12]. 

In the Bjorken kinematic limit, DIS selects the most singular behavior of a correlator of quark fields of the target 
with light-like and causal distance separation z 2 ~ + . With incident photon momentum along the negative 3-axis, 
the kinematics selects z + ~ z±_ ~ leaving z~ as the finite distance conjugate to quark momentum component xP + , 
where x — Q 2 /2P-q is the Bjorken variable, q 2 = —Q 2 is the spacelike virtuality of the photon, and P is the target 
momentum. To leading order in the operator product expansion, the target structure functions are proportional to the 
probability distribution function q/(x). This in turn is given, in manifestly invariant form, by the correlator [45, 46, 47] 

q f (x) = ±f dXe-^{K{P)\%(Xn)Y-n Wf (Q)\K{P)) c . (5) 




X 



FIGURE 1. Left Panel: Pion valence quark distribution evolved to (5.2 GeV) 2 . Solid line is the full DSE-BSA calculation [32]; 
dot-dashed line is the semi-phenomenological DSE-based calculation of Hecht et al. [33]; experimental data points are from [28] 
at scale (4.05 GeV) 2 ; the dashed line is the recent NLO re-analysis of the experimental data [34]. Right Panel: The ratio of u-quark 
distributions in the kaon and pion. The solid line is our preliminary result from DSE-BSE calculations [35, 32, 3 1] ; the experimental 
data is from [36, 26]. 



where / is a flavor label. In the frame where the target has infinite momentum, q/(x) is the quantum mechanical 
probability that a single f-parton has momentum fraction x^k-n/P-n [47]. In the above nt 1 , and partner /?^, are 
light-like vectors satisfying n 2 = p 2 = and n ■ p = 2 and they form a convenient basis for the longitudinal sector of 
4-vectors. The dominant component of q is parallel to n. One can always choose a frame such that k-n — k + and 
k- p = ]C . Note that qf(x) = —qf(—x), and that the valence quark amplitude is q v Ax) = q/(x) — qf(x). It follows from 

Eq. (5) that J Q l dxq v j(x) = (n(P)\Jf(0) -n\n(P)) /2P-n = F K (0) = 1. Approximate treatments should at least preserve 
vector current conservation to automatically obtain the correct normalization for valence quark number. 

In ladder-rainbow truncation, which treats only the valence qq structure of the pion, Eq. (5) yields the explicit form 

rA 

q}(x) = -j ti cd [r % (p,P)S{p)r^(p;x)S(p)r^p,P)S(p-P)] , (6) 

where T^(p;x) is a generalization of the dressed vertex that a zero momentum photon has with a quark. It satisfies 
the usual inhomogeneous BSE integral equation (here with a LR kernel) except that the inhomogeneous term is 
y-nS(p-n-xP-n). This selection of LR dynamics exactly parallels the symmetry -preserving dynamics of the 
corresponding treatment of the pion charge form factor at q 2 = wherein the vector current is conserved by use 
of ladder dynamics at all three vertices and rainbow dynamics for all 3 quark propagators [11, 12]. Here the number of 
valence M-quarks (or d) in the pion is automatically unity since the structure of Eq. (6), along with the canonical 
normalization of the qq BS amplitude r n (p,P), ensures J dxq x u (x) = 1 because / dxT^"\p;x) gives the Ward 
Identity vertex. 

Eq. (6) is in Minkowski metric so as to satisfy the constraint on p ■ n, but LR dynamical information on the various 
non-perturbative elements such as S(p) and F K (p,P) is available only in Euclidean metric [10]. Since q/(x) is obtained 
from the hadron tensor W^ v which in turn can be formulated from the discontinuity T^ v (e) — T^ v (—e), we observe 
that all enclosed singularities from the difference of Wick rotations cancel except for the cut that defines the object 
of interest. With use of numerical solutions for dressed propagators and BS amplitudes that give an accurate account 
of light quark hadrons, our DIS calculations significantly extend the exploratory study made in Ref. [33]. That work 
employed phenomenological parameterizations of these elements. 

In Fig. 1 we display our DSE result for the valence M-quark distribution evolved to Q 2 = (5.2 GeV) 2 in comparison 
with nN Drell-Yan data [28] with a scale quoted as Q 2 = (4.05 GeV) 2 . We also compare with a recent NLO reanalysis 
of the data at scale Q 2 = (5.2 GeV) 2 . The distribution at the model scale Qq is evolved higher by leading order DGLAR 
The model scale is found to be Qq = 0.57 GeV by matching the x" moments for n = 1, 2, 3 to the experimental values 
given independently at (2 GeV) 2 [48]. Our momentum sum rule result f dxx (u a +d a ) = 0.74 (a = 7t),0.76 (a = K) 



at go clearly shows the implicit inclusion of gluons as a dynamical entity in a covariant approach. The large x behavior 
in the form u(x) <~ (1 — x) a has long provided a contrast between previous parameterizations of experiment (a <~ 1.5) 
and QCD-based approaches (a = 2 + y) where /represents anomalous dimensional behavior that increases with scale. 
At our model scale, a = 2.08. Recently, the experimental data [28] has been reanalyzed [49] by an accounting for soft 
gluon processes, and the resulting experimental u(x) at Q 2 = (4 GeV) 2 , forx > 0.4 is essentially identical to the Hecht 
et al. DSE calculation shown in Fig. 1. This eliminates the discrepancy between QCD and experiment here. 

The ratio uk/u k measures the dynamical effect of the local environment. In the kaon, the w-quark is partnered 
with a significantly heavier partner than in the pion and this shifts the probability to relatively lower x in the kaon. 
Our preliminary DSE model calculation [35, 32, 31] is shown in Fig. 1 along with available Drell Yan data [36, 26]. 
Here we include only the leading two invariants of the pion BS amplitude, E(k;P) and F(k;P), where k is qq relative 
momentum. For both amplitudes only the zeroth Chebychev moment in k ■ P is employed . This variable does not occur 
in static quantum mechanics, nor in the Nambu-Jona-Lasinio point-coupling field theory model [38] which also cannot 
accommodate the k 2 dependence. We do not make such a point meson approximation here; the k 2 dependence comes 
from the BSE solutions. Nevertheless, the qualitative features of the ratio uk/u^ in Fig. 1 are adequately reproduced 
by a generalized Nambu-Jona-Lasinio model [31]. 

An explicit analysis of the exponent a in the large-x behavior q(x) ~ (1 — x) a can be made in the limit where the 
propagators in Eq. (6) have the dynamical quark self-energy replaced by a constant (constituent) mass and the vertex 
is taken to be y- n8(p ■ n — xP ■ n). These limits maintain consistency with the Ward Identity mentioned above. We 
also truncate the pion Bethe-Salpeter vertex to just Y$E K {k 2 ) and take E K {k 2 ) = N K /(k 2 — A 2 ) where the relative qq 
momentum is k = p — P/2. For x > 1/2 the pole in the spectator quark propagator is the only one in the upper half 
plane and the p~ integral may be easily evaluated by the residue theorem to yield 

qix) =wL d ^ [x(l-x)P 2 -p\- m 2 ] 2 ] ' (?) 

where now k 2 is evaluated at the p~ pole. It is convenient to change the integration variable to jj. = —p 2 where the 
latter is the value at the p~ pole. Then 

4N c 7t 2 r , xP 2 + ll+m 2 



(8) 



where jJ, m (x) — a/ (I —x)—xm 2 , with a = xm 2 . This divergence of the lower limit for large x renders the result totally 
dependent upon the ultraviolet behavior of the propagators and bound state amplitudes. 

The leading large x behavior of q v (x) ~ (1 — x) a comes from the leading divergent term of each factor of the 
integrand. The integral can be expressed 



^= N L d "[^w^x +d+ ^-" 



(9) 



where the physical case of Bethe-Salpeter bound state amplitudes E 2 determined by one gluon exchange has n = 2. 

The quantities a,b,c,d depend on the mass-dimensioned quantities and scales in the system and have only a 
nonsingular dependence upon x. Their details are not important. A change of variable to p. = ( 1 — x)fi /a shows that the 
integral scales as [(1 — x)/a] n when a/{\ —x) is greater than any physical mass scale in the system. Note that it is the 
pair of bound state vertices that totally determine the large x exponent: if the argument of the bound state amplitude E n 
was not singular at large x, the combined scaling effect of the propagators would vanish giving a = 0. The presence of 
true bound state vertices provides a natural regularization of the integral and allows us to ignore explicit mathematical 
regularization at a higher scale. In a renormalizable field theory such as QCD, there is an overall regularization mass 
scale that can be taken to be much larger than any scale at which physical observable questions are to be asked of the 
formalism. In the above case the renormalization scale is to be much larger than a /(l —x) for values of x where an 
asymptotic exponent is to be quoted. 

If one were to consider the limit of a Nambu-Jona-Lasinio model, the bound state vertices become constants and the 
model integral for q v (x) must be regularized, and the result, including the high x exponent becomes dependent on reg- 
ularization scheme [31]. However, independent of scheme, any thus obtained exponent would be unconstrained 
by the physical properties of the system. 
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